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^ ■ 1 INTRODUCTION 



ABSTRACT 

We construct two linear filtering techniques based on weak gravitational lensing to con- 
strain the inner slope a of the density profile of dark-matter haloes. Both methods combine all 
available information into an estimate of this single number. Under idealised assumptions, a 
is constrained to ^ 15% if the halo concentration c is known, and to 30% if not. We argue 
that the inevitable degeneracy between density -profile slope and halo concentration cannot be 
lifted under realistic conditions, and show by means of Fisher-matrix methods which linear 
combination of a and c is best constrained by our filtering of lensing data. This defines a new 
parameter, called Pi, which is then constrained to ^ 15% for a single massive halo. If the 
signals of many haloes can be stacked, their density profiles should thus be well constrained 
by the linear filters proposed here with the advantage, in contrast with strong lensing analysis, 
to be insensitive to the cluster substructures. 

Key words: Cosmology: theory Cosmology: dark matter Physical data and processes: gravi- 
tational lensing Methods: analytical 



Numerical simulations of non-linear structure formation in a broad 
class of cosmological models, even with different types of power 
spectra for the dark-matter (DM hereafter) density fluctuations, re- 
veal a typical shape for the density profile of DM haloes. As far 
as the numerical resolution allows this statement, the density pro- 
file begins with at least a mild singularity in the core, then falls off 
with a relatively flat slope out to a characteristic radius where it 
gently steepens towards an asymptotic behaviour oc r^'^ far away 
from the core. Do real haloes behave in the same way as theory 
predicts? 

Gravitational lensing should in principle be able to give the 
cleanest answer to this question. Density profiles in galaxy-sized 
objects are expected to be modified on small scales by baryonic 
physics, where they are likely to approach the isothermal den- 



The question is important because it aims at a central predic- 
tion of non-linear cosmological structure formation. Answering it is 
complicated by the angular resolution limit of ^ 20" of weak grav- 
itational lensing, set by the number density of background galaxies, 
and by the high non-linearity of strong gravitational lensing. In fact, 
claims that strong gravitational lensing, when combine d with stel- 



sity slope cx 



instead of the generic DM behaviour. On the 



mass scale of galaxy groups or clusters, however, baryonic physics 
should be constrained to the innermost region, leaving the DM den- 
sity profile almost intact. Galaxy-gala xy lensing seems to sho w 
tentative evidence for this expectation jMandelbaum et al]l2006l) : 
while the shear profile around low-mass haloes is consistent with 
an isothermal density profile, it seems to flatten towards the theo- 
retical expectation for DM haloes around high-mass haloes. 



lar dy namics, requires flat halo c ores have been made i Sand et afl 
12004) and doubted. In particular jMeneghetti et alj|2007l) showed 
how the measurement of the inner slope can be systematically un- 
derestimated if halo's substructure are not taken into account. A 
weak lensing analysis, even if observationally more challenging, 
has the advantage to be almost insensitive to cluster's substructures 
because of the instrinsic nature of the signal. 

Previous studies based on weak lensing have followed an ap- 
proach where a shear profile was first measured and then fit to 
the shear profile expected from certain three-dimensional density 
profiles, thus indirectly constraining the density-profile models. 
Given the sparseness of lensing information near the core of galaxy 
groups and clusters, we develop a different approach here. Instead 
of constraining the shear profile, we only wish to derive a single 
number from the shear data, namely the slope a of the density pro- 
file within the characteristic radius, assuming that the asymptotic 
outer slope is —3. 
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Figure 1. Shear profile of a generalized NFW halo for three different values 
of a. The case a = 1 corresponds to the usual NFW profile. Note the strong 
dependence on the inner slope. 



We pursue this approach with two linear filtering techniques. 
One of them is specifically constructed below to return a as its 
only result. It is thus made to combine all available information 
into its estimate and should thus optimise the significance of the 
measurement. The other varies the inner slope of the density profile 
until it finds the maximum signal-to-noise ratio in a given sample 
of haloes. 

We proceed as follows. In Sect. 2, we introduce the weak- 
lensing properties of the generalised NFW density profile. We de- 
velop our linear filters in Sect. 3. There, we also discuss the de- 
generacy between the central density slope and the concentration 
parameter of the halo and use the Fisher matrix to find a linear pa- 
rameter combination which is best constrained by shear measure- 
ments. The sensitivity of our method and its limitations are shown 
and discussed in Sects. 4 and 5, respectively. Section 6 presents our 
conclusions. 



2 GENERALIZED NFW PROFILE 

In the last decade a large effort has been devoted to predict the 
density profile of DM haloes in KCDM cosmologies. Thanks 
to increasing resolution of numerical simulations there is now a 

' It r * n ! 

common agreement (Navarro et al. 1997; Moore et al. 1998), that 
DM density profiles can be accurately described by the generalized 
NFW profile 

Ps 



^^'■^ (r/r,)«(l + r/r,)3— 
with the scale radius r^, the inner slope a and the scale density 
200(3 -a)(r2oo/rs) 



Ps — pcrit{z) 



3j-F2(3 — a, 3 — a, 4 — a, — r2oo/rs) 



(1) 



(2) 



where \F2{a, b, c, z) is a hypergeometric function and r2oo is the 
radius enclosing 200 times the critical density of the universe pcrit- 
The halo concentration is defined as 



C200 = 



''200 



(3) 



Following iKeeton & Madaj ( 1200 ih we interpret the scale ra- 
dius as the radius where the density profile reaches slope —2, 



i.e. dlnp/dlnr 



and thus 



-2. For the profile of equation[T] 

r-2 = rs{2- a) 



C-2 



'"200 

r-2 



-C200. 



(4) 



(5) 



For a = 1 this f ormula reduces to the standard NFW case 
dNavarro et al. I ll997h . The profile is fully characterized when the 
mass, the redshift, the concentration and the inner slope of the halo 
are specified. However not all of these parameters are independent. 
Numerical simulations show that it is possible to define fitting for- 
mulae relating the concentration with the mass and the redshift of 
the halo. In the f ollowing we will use the prescription proposed by 
lEke et al.l ( l200lh for this purpose. They also found that for a fixed 
value of mass and redshift, the concentration approximately follows 
a log-normal distribution 



p{c)dc — 



exp 



(Inc — Inc) 



2cr2 



din c 



(6) 



where is the l-a dev iation of A (Inc) ~ 0.2 dNavarro et al.l 
ll997l;lBullock et alj200ll) . 



2.1 Weak lensing properties 

This Section summarises the basic weak-lensing concepts that 
will be used later. Fo r a c omplete overview we refer to 
iBartelmann & SchneideJ ( 1200 ll) . An isolated lens with surface 
mass density E(0) has the lensing potential 



*(0) = 



4.G DiDs 

c2 Dls 



<re E(6» )ln|0-0 



(7) 



where G is the gravitational constant, c is the speed of light, and 
Di^s,is are the angular diameter distances between the observer and 
the lens, the observer and the source and the lens and the source 
respectively. 

Due to the presence of the lens a light ray is deflected by the 
angle 



C(0) = V*(0). 



(8) 



A source located at the angular position /3 in the sky is seen by 
the observer at an angular position 6 which are related by the lens 
equation 



f3 = o-m- 



(9) 



If the source's angular extent is much smaller than the angular scale 
on which the lens properties change, the lens mapping can be lo- 
cally linearised and the image distortion is given by the following 
Jacobian matrix 



^- do 



where 



9^^(6>) 



1 — K — 71 
—72 



—72 

1 - K + 71 



K(0) = ^ = i(*n + *22) 



(10) 



(11) 



is the convergence, i.e. the surface mass density scaled by the crit- 
ical surface mass density 



Ds 



AtvG DiDis ' 



(12) 
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and 



71 = - *22),72 = *i 



(13) 



are the two components of the complex shear 7 = 71 +172 • The net 
resuU is a distortion and a magnification of the baclcground sources 
due to the lens gravitational field. 

For an axially-symmetric lens, outside critical curves, the dis- 
tortion is tangential to the line connecting the source and the lens 
so that the tangential shear is given by 

7T = -(71 cos(26l) + 72 sin(26l)) = -^(76^'*) (14) 

and the shear modulus can be derived from the convergence, 

\-i\{e) = R{9) - K{e) (15) 

where R{6) is the mean surface mass density inside a circle of ra- 
dius 6 centered on the lens and hi{6) is the convergence at radius 
6. If a 7^ 1 it is not possible to find an analytic expression for 
the shear profile and therefore equation [T5]has to be computed nu- 
merically. We show in figure [T]the shear profile for three values of 
a (0.5,1.0,1.5). The inner shear profile depends sensitively on the 
inner slope a. In the case a = 1 the shear is logarithmically di- 
vergent for small values of e ( lBartelmanrj[T996i) . The divergence 
is more pronounced for steeper profiles, while the shear profile de- 
creases for Q < 1 and tends to converge to a finite value even if it 
is undefined for ^ = 0. 

The shear profile depends on the two parameters of the den- 
sity profile, of which the concentration depends mildly on the halo 
redshift. An additional and stronger dependence on halo and source 
redshifts is introduced through the geometry of the lens system. 

The observable lensing signal is the ellipticity of the back- 
ground galaxies. 



+ 5 



1 + g*e^' 



where 



h/a 



exp(2i0) 



(16) 



(17) 



1 + b/a 

is their intrinsic ellipticity, a and h are the major and minor axes, 
respectively, is the orientation angle, g is the reduced shear 
9 — l7l/(l ~ 9* is its complex conjugate. In the weak- 

lensing regime (7 ^ 1), e° " ~ p ~ 7. In the following, we will 
exclusively use the reduced shear since we want to explore scales 
where the approximation k ^ 1 does not hold. Nonetheless, we 
shall denote it by 7 throughout for simplicity of notation. 



3 METHODS TO CHARACTERIZE THE SHEAR 
PROFILE 

In this section we describe two methods, based on optimal linear 
filters jSanz et al. 2001; Maturi et al. 2005), to estimate the inner 
slope of DM haloes using weak-lensing observations. The advan- 
tage of linear filtering as opposed to standard profile fitting is that 
filters can be constructed such as to minimise noise caused by in- 
tervening structures along the line-of-sight. 

3.1 Optimal linear filtering 

For a generic optimal linear filter, the data D{0) is modelled as the 
sum of the signal to be measured and the noise 

D{e) ^ S{9) + N{9), (18) 



where S{0) = At{G), A is the signal amplitude and t{G) is a 
model for its angular shape. In our application, the signal is the 
lensing shear of the intervening DM halo and the noise is given 
by the intrinsic ellipticity of the background galaxies, their finite 
number and the contamination due to large-scale structures. The 
noise components are assumed to be Gaussian, random with zero 
mean and isotropic since their statistical prope rties are independe nt 
of the position in the sky (for further detail see (' Maturi et alj2 005S). 
We now define a linear filter ^'(0, a, w) which, when convolved 
with the data, yields an estimate for the amplitude of the signal at 
the position G: 



Aest(G) 



which is unbiased 



A 



D{G')-^{e-G',a,w)d^e', 



•^{G,a,'w)T{0,a,w)d^e ~ 1 



and whose variance cr 



27r 



\^{k,a,w)\^PN{k)d^k 



(19) 



(20) 



(21) 



is minimal. The filter ^ satisfying these two conditions minimises 
the Lagrangian L — -\- \b. It reads 



*(fc) 



1 

2^ 



|f(fc, a, w) 
Piv(fc) 



f (fc, a, w) 

PNik) 



(22) 



where w — (c, M, z) and VE* and f are the Fourier transforms of 
the filter and the signal shape, respectively . Note that we have as- 
sumed in the previous derivation that the mean values of the halo 
parameters (to) are well known. This is an idealising assumption 
and we refer to Sect. 14. II for a more detailed discussion. The filter 
depends only on the angular shape of the signal r(fc, a, w) and the 
noise power spectrum Pn- In particular it is most sensitive to those 
spatial frequencies for which the signal t is large and the noise 
power spectrum is small. This filter is optimal in the sense that it 
maximises the signal-to-noise ratio for the a given assumed signal 
shape. 

The left panel of Fig.[2]shows the filter's shape calculated us- 
ing three different values of the inner slope, a = 0.7, 1.0, 1.3. 

3.1.1 Dealing with non-linear signaLs 

The filter described in the previous section can be used to measure 
quantities which appear linearly in Eq.[T8](e.g. the amplitude of the 
shear signal). This is not the case for the inner slope a breaking the 
main assumption on which the linear filter is based on. However, 
if we expand the halo's shear profile around a fiducial value of the 
inner slope, ao. 



j{G, a, w) = 7(0, ao,w) + 
Eg.lISIreads 



dj{G, a, w) 



da 



Aa, 



(23) 



DiG)-'y{G,ao,w) 



dj{0, a, w) 



dc 



Aa + N{e), (24) 



such that Aa appears linearly and the linear filtering scheme can be 
applied. The shear derivative with respect to a plays the role of the 
signal shape, r, and Aa that of the amplitude A to be measured. 
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Figure 2. Optimal linear filter (left panel) and scale-adaptive filter (right panels) shapes calculated for three different values of a. All filters are normalised to 
unity. 




Figure 3. Estimated inner slope otest as a function of the fiducial inner 
slope ao used in the linear filter with the l-cr en'or bars calculated via 
Eq.|2T] The black line shows the real value of the halo's inner slope. 



This allows the definition of the following estimator for the inner 
slope, 



where 



A7(0, Q, Qo, ly)) = 7(0,Q,iu) - 7(0,Qo,iu). 



(25) 



(26) 



The approximation applied in Eq. |24] implies that a''" is a 
good estimator of the inner slope only when ao is close to the real 
value of a. If this is not the case, the value of the inner slope tends 
to be overestimated as we show in Fig. [3] If a single halo is con- 
sidered, the error bars associated to aest are so large that the over- 
estimation can be neglected for a large range of qq. However, if 
several haloes are stacked, the error bars shrink and the overesti- 
mation becomes important. In order to avoid this problem, more 
measurements of the same halo have to be carried out sequentially: 
the first measurement starting with an arbitrary value of ao, and the 
second using the estimate a*^"' found previously as a fiducial value. 



We tested that, for a reasonable guess of the first fiducial value, 2-3 
measurements suffice to recover the correct slope. 



3.2 Scale-adaptive filter 

The linear expansion used in the previous section can be avoided 
by defining a scale-adaptive filter. Such a filter is defined similarly 
as the linear filter from Sect. l3.ll with an additional constraint on 
the amplitude of the signal A^st which must be maximised when 
the adopted inner slope fits the data best. 



da 



The minimisation of L = cr^ + Ai6 + leads to the filter 



*(fe, a) 



1 f(fc,Q) 1 



27r Pjv(fc) A 
with the constants 



2fe + c-(2a + 6) 



d In f (fc, a) 
din a 



211 J PN{k) 



27r J PN{k) dlna 



2tv J Pjv(A:) I dlna 



(27) 



(28) 



(29) 



(30) 



(31) 



A = ac-b. (32) 

Its defining property is thus to maximise the signal-to-noise 
ratio when the correct inner slope is adopted. This implies that the 
inner slope can only be determined indirectly from a sequence of 
measurements of the shear amplitude Aest, searching for that value 
of a that maximises A^st- 

The filter shape is plotted in the right panel of Fig. |2] 
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3.3 Dealing with parameter degeneracies 

The two methods presented in Sects. [TT] and 13.2] assume a clus- 
ter model with known mass, redshift and concentration. In a re- 
alistic situation, we can assume to have sufficiently precise red- 
shifts. Mass estimates would have to be obtained from optical rich- 
ness, kinematics of the cluster galaxies or X-ray scaling relations. 
Then, estimates for the concentration could be derived from the 
mass-concentration relation found in numerical simulations, albeit 
with a considerable scatter. The concentration depends only very 
weakly on the mass, hence uncertainties in the mass estimate do 
not strongly affect the concentration estimate, and thus the mass 
does not need to be precisely known. However, numerical simula- 
tions suggest a log-normal distribution of the concentration around 
its mean with a standard deviation of ~ 0.2, which implies that 
concentration parameters of real clusters can only be very poorly 
guessed. 

Moreover, the inner slope, as the parameter we are aiming to 
measure, is degenerate with the concentration. In fact, it is possi- 
ble to describe a halo with high central density with a large value 
of ct and a smal l value of e or vice versa, and so the problem is 
not well defined jWvithe et al.l200lh . Thus, any attempt at measur- 
ing the profile's inner slope depends critically on the assumed halo 
concentration, which is uncertain in reality. 

To cope with this problem, it is convenient to re-parametrise 
the profile accounting for this model degeneracy, defining new pa- 
rameters which can be more precisely measured. In short, the logic 
behind the procedure described below is as follows. In a realistic 
situation, we have no chance to break the degeneracy between c 
and a. Rather, we can rotate the parameter space such that one of 
its axes becomes parallel to the degeneracy direction and the other 
perpendicular to it. The latter will define a new parameter as a lin- 
ear combination of c and a which observations can constrain best. 
Comparisons with theory should then be performed on the basis of 
this parameter rather than through c and a separately. 

This is achieved by a Fisher-matrix analysis. The Fisher ma- 
trix is 



(33) 



where L is the logarithm of the likelihood function and tt = (a, c) 
are the free model parameters. In case of a Gaussian probability 
distribution, the Fisher matrix can be written as 



F,j = -Tr[A^Aj + C'HUj] (34) 

where C is the covariance matrix, Ai — C~^C,i and Mij = 
and ^ is the assumed model. Since C does not depend 
on the inner slope and on the concentration, the first term in Eg. 1341 
vanishes. We evaluate the Fisher matrix at a fiducial point (qq , co ) . 
In particular, we assume uq — 1 and we calculate cq using the pre- 
scription by (Eke et al. 2001). We truncate the shear profile at an 
inner radius rmin = 1/ ^/ngai, which is the minimum achievable 
resolution for a given number density Ugai of background galaxies 
and at an outer radius Tout = '"200 • Once Tout > rs the Fisher 
matrix depends negligibly on rout since the derivative of the shear 
profile with respect to alpha is zero and the derivative with respect 
to the concentration is very small. 

The eigenvectors {vi,V2) and (u3,W4), of the Fisher matrix, 
determining the directions of largest and smallest degeneracy be- 
tween the parameters a and c, define a rotation of the parameter 



space and thus two new parameters 

Pi = Via + V2C 
P2 = vsa + V4C 



(35) 
(36) 



which are linear combinations of a and c. The two new parameters 
are those which can be constrained best and worst, respectively, 
given the model adopted in the Fisher-matrix estimate. 

If the linear filter is used to measure the inner slope, the model 

fj, is 

H = aest{a,c) = j ['y{a,c) - 'y{ao,co)\^{ao,co)(fx (37) 
d/J, _ f d'y{x, tt) 



and thus 



-*(a;, -Kojd X. 



(38) 



The covariance matrix reduces in this case to the variance of the 
measurement obtained from Eg. 1211 

Note that the Fisher matrix defined above is singular, i.e. its 
determinant vanishes. The errors on the new parameters are given 
by l/^fXi, where Ai are the eigenvalues of the Fisher matrix. Since 
one of them is vanishing the error on one parameter (taken to be P2) 
is infinite. This means that the likelihood region in the plane (q, c) 
is an ellipse infinitely elongated in the degeneracy direction. This 
is because there is more than one way of fitting a single data set 
(Aa) by varying the two parameters. In the right panel of Fig.|4]we 
show the result for a halo of A/ = 5 x lO^^A/© /h at redshift z = 
0.3 with concentration c = 4.4. The corresponding eigenvector 
components are «i = 114 = 0.95 and 112 = —vs = 0.30. 

When the scale-adaptive filter is used, the measurable quantity 
is the shear amplitude 



A{a, c) 



D{aH,CH;e)'^{a,co;e)d' 



(39) 



and the value of the inner slope (cvest) is then estimated looking 
for the value of a maximising the amplitude. It is clear that it de- 
pends only on the halo's concentration co assumed in the filter. To 
find the degeneracy direction between the inner slope and the con- 
centration in this case, we analyse the relation between a^st and 
c around a fiducial point in the (q, c) plane. The result is shown 
in the left panel of Fig. |4] for the same halo as considered before. 
Here, too, we define two new parameters Pi = 0.97a + 0.22c 
and P2 = —0.22a + 0.97c. In this case, the error cannot be cal- 
culated analytically since the measurement is indirect. Instead, we 
have performed a Monte-Carlo simulation (see Sect.|4j. 

Since the shapes of the filters are different, so are the degener- 
acy directions we find. 

The probability distributions of Pi and P2 can be found con- 
volving the probability distributions of the concentration and the 
inner slope. Using the degeneracy direction found for the linear fil- 
ter and assumi ng a log-normal di stribution for the concentration 
with ac = 0.2 teullock et al.l200lh and a Gaussian distribution for 
the inner slope with aa = 0.15 jPiemand e t al. 2004), we find that 
both probability distributions of Pi and P2 can be approximated as 
log-normal distributions with standard deviations crpj =0.29 and 
(JP2 = 0.33 respectively, as shown in Fig.[5] 



4 METHOD UNCERTAINTIES 

Here, we discuss in detail possible error sources affecting the 
measurement of the inner slope using the methods described in 



6 



Viola et al. 



Linear filter 



Scale-adaptive filter 
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Figure 4. l-cr and 2-a likelihood regions in the plane (a, c) computed using the linear filter (left panel) and the scale-adaptive filter (right panel) for a halo of 
M = 5 X lO^Mg/h at ^ = 0.3. The fiducial value is (1.0, 4.4). 



Distribution function 



Effective number counts 




Figure 5. Probability distributions for a, c, Pi and P2. A normal distri- 
bution with CTq = 0.15 for the inner slope jPiemand et a l and a 
log-n ormal distribution with = 0-2 for the concentration (Bullock et al. 
I2OOII) are assumed. The probability distribution for Pi is given by the con- 
volution of the probability distribution of a and c, while the probability 
distribution for P2 is given by their cross correlation. 



Sects. 13.11 and 13.21 We will show the error calculation for a halo 
of M = 5 X IQ^'^Mq/K z = 0.3, c = 4.4. 

The statistical uncertainties arising from the data noise com- 
ponent A'^ are given by the intrinsic ellipticity of the background 
galaxies, their finite number and from the contamination due to 
the intervening large-scale structures. The filters we have defined 
minimise these uncertainties. They are quantified by Eq.[2T|for the 
linear filter and by a Monte-Carlo analysis for the scale-adaptive 
filter since in this case ct is measured indirectly by estimating the 
location of the maximum in the estimated signal, and an analytical 
computation of its variance is impossible. 

The Monte-Carlo analysis has been performed generating 
1000 realisations of a shear catalogue using randomly distributed 
background galaxies with a density Ugai ~ SO/arcmin^ , placed 
at redshift Zs — 1.0, on a 0.01 degree field. The halo has been 
placed in the field center. The noise due to the intrinsic galaxy el- 
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Figure 6. Magnification bias expected for an halo of 5 X lO^'^M©//! at 
redshfit z = 0.3 tensing galaxies aX z = 1.0. 7 is the exponent of the 
power low in equation |40l 



lipticities ((Te = 0.3) and the lensing effect due to the intervening 
large-scale structure have been added. The latter noise is calculated 
assuming that the large-scale structure can be described by a Gaus- 
sian random field with a power spectrum determined by the linear 
theory of structure growth. 

We assume in our analysis that the magnification bias can be 
neglected, allowing us to leave the effective number ngai of avail- 
able galaxies unchanged. This is justified only if the slope 7 of the 
flux distribution of faint galaxies 

no(> S) = aS"^ (40) 

is unity as discussed by iBartelmann & Schneide3 ( 1200 ih . The ef- 
fective number of galaxies n^f f scales with 7 as 

= (41) 

ngai(> k>) 

where /i is the magnification. Specifically, n^// is lowered by at 
most 40% compared to Ugai near r = 0.2rs if 7 is 0.5, as shown in 
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Figj6] For galaxies in the Hubble Ultra Deep Field l lBeckwith et alj 
l2006h we estimate 7 ~ 0.8 causing a magnfication bias of around 
10%. 

For each realisation we use Eq.[T9]to estimate the shear ampli- 
tude in the position corresponding to the halo's center using filters 
initialised with an inner slope in the range [0.6 — 1.4]. The esti- 
mated inner slope value is then defined as the value of a giving the 
maximum value of the shear amplitude. We finally calculate their 
distribution and the dispersion around the mean value (the results 
are summarised in the fourth column of Table 1). 

We find that the standard deviation associated with the inner 
slope, measured by the scale-adaptive filter, is 0.19. The analytical 
calculation done for the linear filter gives a value of 0. 14. 

The same calculation has been done considering haloes of dif- 
ferent masses and at different redshifts. As shown in Fig. |9l the 
standard deviation increases with respect to the redshift and de- 
creases when the mass is increasing. In particular for a halo placed 
at intermediate redshift between the background sources and the 
observer, the standard deviation varies in the range [0.2 — 1.0] for 
amass range [10^^ - 5 x lO"]. 

The preceding calculations show that errors on the inner slope 
due to intrinsic ellipticities of background galaxies and due to con- 
tamination by large-scale structures are large when computed for a 
single halo. However stacking a large number of haloes (10-100), 
it is possible to measure an average value of a with a few percent 
accuracy. 

A more accurate error evaluation has to consider also the 
scatter around the fiducial value of the halo's mass, redshift and 
concentration used in the filter definition. For both methods, we 
perform a Monte-Carlo simulation, following the procedure de- 
scribed above, assuming a Gaussian distribution for the halo mass 
(o"A/ = 1.5 X 10^*) and redshift (ctz = 0.03) and a log-normal 
distribution for the concentration (ctc = 0.2) following numerical 
simulations (citation). The result is shown in Fig.|7] 

One critical point that we have avoided so far concerns the 
choice of the fiducial values for the haloes parameter. We discuss 
this point in the following section. 



4.1 Model sensitivity 

Defining the filter requires the specification of a model. The estima- 
tor (Eq.[79l) we defined for the inner slope is unbiased only if the 
model is correct. We investigate here what happens if the filter is 
defined using a generalised NFW profile with wrong fiducial values 
of mass, redshift and concentration. We study in particular the case 
in which the fiducial redshift used in the filter differs from the real 
redshift by about 10 %, the mass by about 30 %, and the concen- 
tration by about 20 %. We show the results in the first three panels 
of Figs. [T2] and [T3] (blue lines) for the linear and the scale-adaptive 
filter, respectively. 

As expected, the inner-slope estimate is biased. This reflects 
the degeneracy between the parameters, in particular between the 
scale radius r_2 = r2W) (M, z)/c-2 and the inner slope. The scale 
radius depends only slightly on the halo mass and redshift, while it 
is strongly affected by a variation in the concentration. 

This bias has to be compared with the statistical errors associ- 
ated with the measurement in order to assess whether uncertainties 
in the fiducial halo parameters are important or not. If a single halo 
is considered, a wrong assumption on the concentration (the most 
critical parameter) introduces a bias that is on the same order as the 
statistical error. However, if several haloes are stacked (we show in 



Fig.[T2]results after stacking 10 and 100 haloes), the bias is a factor 
of 10 larger than the statistical uncertainty. 

In Sect. 13.31 we discussed how it is possible to deal with de- 
generacies between inner slope and concentration, defining two 
new parameters (Pi, P2), linear combinations of c and a, which 
are respectively the best and the worst constrained parameters given 
our model. The measurement of the new parameter Pi is almost un- 
affected by the choice of the other parameter P2 as we show in the 
right panel of Fig.[T2]and[T3] while the effect of a wrong assump- 
tion of halo mass and redshift produces a similar bias. We recall 
that these latter quantities can be measured by means of other ob- 
servables, as discussed before. 

Once the model had been re-parametrised in term of Pi and 
P2, we estimated the error on Pi using a Monte-Carlo simulation 
in the same way we have done before for a. The result is shown in 
Fig.i 



5 POTENTIAL PROBLEMS 

We now want to point out the conditions under which the two meth- 
ods described can be successfully applied. 

First of all, the reduced shear must be measurable at relatively 
small angular scales (smaller than the scale radius of the halo) 
where the density profile is sensitive to a change of the inner slope. 

Towards the halo's centre, the image distortion becomes non- 
linear such that the galaxy ellipticities are no longer an unbiased es- 
timator of the shear. We quantify the expected deviation by a simple 
test: We use the deflection-angle map of an NFW halo to lens a cir- 
cular source (for which we assumed a Sersic profile with n = 1.5 
and r = 0.35arcsec) moving radially towards the halo centre. We 
measure the ellipticity of its image (using quadrupole moments) as 
a function of cluster-centric distance and compare it to the true re- 
duced shear Figure [TO] shows the result for three different haloes 
{M = 10", 5 X 10", 10^^ A/0 //i). The conclusion is that up to 
r = 0.2rs the measured ellipticity of galaxies is still an unbiased 
estimator of the (reduced) shear while at smaller scales the con- 
tribution from higher order terms start to be dominant. Therefore, 
r ~ 0.2rs should be taken as the minimum radius where the mea- 
sured ellipticity can still be considered to faithfully represent the 
reduced shear. 

However measuring shear at these scales can be tricky even 
with a high background galaxies density due to the possible di- 
lution of the shear signal caused by cluster galaxies. In order to 
avoid this problem, accurate colour-magnitude information should 
be available so that it is possible to well se parate cluster members 
from non-members l lBroadhurst et alj|2005l) . 

We showed in the previous section that the error associated to 
the measurement of the inner slope is high when computed for a 
single halo. Thus several haloes need to be stacked together The 
number of haloes to be stacked depends strongly on the minimal 
radius where the shear can be detected, and on the number of back- 
ground galaxies. In Fig.[TT] we plot the relative error on the mea- 
surement of a as a function of these two parameters for a halo of 
M = 5 X IO^^A/q/Zi at redshift z = 0.3. Assuming 30 galaxies 
per square arc minute, the number of haloes to be stacked to reach 
an accuracy of a few percent on the inner slope is between 10 and 
100 going from Tmin ~ 0.2rs to rmm = O.Sr^ . We emphasise that 
the stacking procedure can be affected by a wrong determination of 
the cluster centre that causes a circularisation of the average cluster 



profi le in its central part |Kathinka Dalland Evans & B ridle 2008). 
iMeneghetti et al.l j2007h showed how the determination of the 
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Figure 7. Normalised distributions of the value of the inner slope computed using the linear filter (left panel) and scale-adaptive filter (right panel). A 
Gaussian distribution has been assumed for the mass and the redshift, while a log-normal distribution has been adopted for the concentration. 
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Figure 8. As Fig.]?] but using the new parameters Pi and P2 instead of a and c. 



Table 1. Statistical eiTors in the parameters measurement for a halo of 5 X 10^^ /h at redshift z = 0.3. In the first column we indicate the used filter, in the 
second column the parameter we constrain and in the third column its fiducial value. In the fijurth column are shown the expected errors assuming randomly 
distributed background galaxies with intrinsic ellipticity CTe = 0.3 and random noise due to the large-scale structures. The errors presented in the fifth column 
take also into account Gaussian errors in the halo mass and redshift with standard deviations (Tm = 1.5 x 10^* and Uz = 0.03, respectively, and a log-normal 
distribution for the concentration with standard deviation ac = 0.2. When Pi is estimated the probability distribution of P2 is calculated from the probability 
distribution of the concentration assuming a Gaussian probability distribution for the inner slope with CTq = 0.15. In the sixth column we show the percentage 
error on the parameter estimation. 



Filter 


Parameter 


Fiducial value 


cr (stat.) 


(T (stat.-l-model) 


Percentage error 


SAF 
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1.00 


0.19 
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PI 


1.95 


0.21 


0.29 


0.15 


LF 
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1.00 


0.14 
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LF 


Pi 


2.31 
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0.13 
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Figure 9. Standard deviation for o as a function of the halo's redshift and 
mass. 
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Figure 10. Comparison between the theoretical reduced tangential shear 
(thick lines) and the shear estimated from galaxy ellipticities (thin lines) for 
three different masses. 



inner slope can be biased if the triaxiality structure of the haloes 
are not taken properly into account. However if many haloes are 
stacked together a direct comparison with the projected DM aver- 
age profile found using stacked simulated clusters can be consis- 
tently done. 

Moreover the effect of the baryons in shaping the density pro- 
file at this scale is not negligible. We plan to attack this problem 
using numerical simulation in order to study the effect of stacking 
and the presence of the baryons on our results. 



6 CONCLUSION 

Starting from the question how the central density profiles of group 
or cluster-sized, DM haloes can best be constrained and compared 
to observations, we have developed two methods based on linear 
filtering of gravitational-shear data that aim at returning a single 
number, i.e. an estimate of the inner slope a of density profile. One 
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Figure 11. Standard deviation of the inner slope as a function of the 
background-galaxy number density and of the minimal radius where the 
shear can be detected. The calculation has been done for a halo of M = 
5 X IO^^'Mq/Zi at redshift z = 0.3 with a = 1.0 using the linear filter. 



filter is constructed to directly return this number, the other searches 
for the maximum of the signal-to-noise ratio as a function of a. Our 
results are as follows: 

• When applied to a single halo of 5 x 10^^ Mq /h near z — 0.3, 
the inner slope of the density profile can be estimated with a I-ct 
accuracy of 14 % with the linear filter and 19 % with the scale- 
adaptive filter, provided the halo concentration is known. Even 
though this situation is unrealistically idealised, it is promising be- 
cause it is based on a single halo only. 

• Taking the considerable uncertainty in halo concentrations 
into account increases the l-cr error to between 25 ... 30 %. 

• In reality, the halo concentration is at best roughly known. 
Based on real data, there is an almost perfect degeneracy between 
a and the halo-concentration parameter c: if c is assumed to be too 
large, a will be underestimated and vice versa. Based on lensing 
data only, this degeneracy cannot be lifted. 

• To address this problem, we search for that combination of the 
parameters a and c that can best be constrained by observations. 
We set up the Fisher matrix, rotate the two-dimensional parameter 
space to diagonalise it and identify its smaller eigenvalue as that 
best-constrained parameter, called Pi. We find Pi = 0.95q+0.30c 
for the linear filter and Pi = 0.97a + 0.22c for the scale-adaptive 
filter. 

• This parameter Pi is now constrained with a l-a relative ac- 
curacy of ~ 14 % both with the linear and the scale-adaptive filters 
and the measurement is almost insensitive to the value of the other 
parameter P2. 

While these results seem highly promising, in particular when 
applications to cluster samples rather than individual clusters are 
envisaged, we consider our study as a first step. While we have 
taken into account that image ellipticities measure the reduced 
gravitational shear rather than the shear itself, measuring the re- 
duced shear near the centres of galaxy groups or clusters is severely 
hampered by the cluster galaxies themselves. It thus appears nec- 
essary to stack the signal from several or many clusters to arrive at 
a reliable estimate for a. Then, clusters with different masses, red- 
shifts and concentration parameters will inevitably be combined. 
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Figure 12. Left panels: Estimated inner slope of the halo (aest) as a function of the fiducial inner slope used in the filter (oq) with the 1 — cr 
error calculated for 10 (orange) or 100 (green line) haloes using the Monte-Carlo simulations described in Sect. |4] The black line shows the 
real value of the halo's inner slope. The first panel shows the bias caused by a fiducial concentration 20 % larger or smaller than the real 
concentration. The second panel shows the bias induced by a 50 % difference between the fiducial and the real halo's mass, while the third 
panel shows the bias caused by a difference of 10 % between the fiducial and the real halo's redshift. Right panels: As the left panels, but for 
the new pair of parameters P\ and P2 . 



with the tendency to blur the signal. However, the resuhs derived 
and presented above indicate that the principle of our approach is 
promising, which consists in combining all available information 
into a single number, which is thus well constrained. Further stud- 
ies are required to address the issue of stacking data in this context. 
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